Can dileptons be observed in heavy ion collisions at RHIC? 
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J>V Abstract 
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Both dilepton and charm production at RHIC are considered to be important 
signatures for Quark-Gluon Plasma production. Recently it was argued by S.Gavin, 
P.L.McGaugheyP.V. Ruuskanen and R.Vogt ||] that the background from semilep- 
tonic correlated charm decays is so large that it makes dilepton measurements vir- 
tually impossible. We show that this conclusion in fact is in fact reversed if the 
energy loss due to secondary interaction of charmed quarks is included. 

O ' 

Dileptons produced in highly excited hadronic matter provide valuable information 
about the hottest and the most dense stages of nuclear collisions. Together with photons, 
they are the so called penetrating probes which suffer very little secondary interac- 
tion. Consequently, dilepton measurements has attracted great deal of attention, both 
of theorists and experimentalists. Referring specifically to highest energies, let us re- 
mind that one of the major RHIC detectors, PHENIX, and ALICE at LHC plan dilepton 
measurements, both with electron and muon pairs. 

Another potential QGP signature is thermal production of new quark flavors, especially 
of charm [@, |J [5|] [] However, the ordinary partonic production of charm at the first impact 
is large, and whether it dominates the secondary (and thermal) charm production remains 
unclear. Nevertheless, charm signal will also be experimentally addressed at RHIC, by 
STAR and PHENIX collaborations. 

Since charmed hadrons have substantial semileptonic decays, their simultaneous decay 
create a background for dilepton measurements. This issue was addressed recently 
by S.Gavin, P.L.McGaughey,P.V. Ruuskanen and R.Vogt (below GMRV) in a detailed 
paper 0. Their conclusions are summarized in Fig]l|, and they basically imply that the 
background from leptonic decays of charmed (and even b) quarks is so large, that imple- 
mentation of the dilepton measurements is virtually impossible in the whole kinematic 
domain. 

In this short paper we question those pessimistic conclusions, and suggest that it 
should be reversed. We show that a very important effect is missing from the GMRV 

1 The expected highest temperatures at RHIC |3[ |) are Ti = 400 — 500Af eV, so the mean energy 
per parton w 3T is comparable to charm quark mass. Note also that the mass of the strange quark is 
not large enough to suppress its production in the hadronic phase. 
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Figure 1 




Figure 1: Contributions of different dilepton production mechanisms according to GMRV 
for central Au Au collisions at RHIC (a) and LHC (b). The curves correspond to: a cor- 
related charm decays (dash-dots), b-quark decays (dash-double dots), Drell-Yan process 
(dashed curve), thermal dileptons (solid curve) and decays of thermally produced charm 
(dots). 

analysis: unlike dileptons, the charmed quarks are not "penetrating probes", and their 
spectra should be very different in pp and heavy ion collisions. Like any other quarks 
(and gluons), charmed ones are also subject to energy losses due to multiple secondary 
interactions in dense matter produced in the collisions. As we will show below, they are 
mostly stopped in matter. 

Our first point is a purely geometrical observation. Consider collision of two heavy 
nuclei, and imagine that in it a cc pair is produced^. The charmed quarks have to pass 
certain distances di, d 2 on their way out, and we point out that it is very improbable that 
the sum di + d 2 is small because the quarks are mostly produced back-to-back. 

If nuclei are approximated as spheres with a well-defined surface and radius R, one 
can easily quantify the relevant distributions. A distribution over a single quark path 

2 The commonly used terminology a correlated and an uncorrelated charm decay. The former is a 
simultaneous decay into l + and l~ from a cc pair produced in one parton collisions, while the latter 
comes from the charm quarks produced independently. In this note we concentrate on the correlated 
background only because the uncorrelated background can be statistically subtracted in a standard way. 
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Figure 2: The histogram shows distribution of the transverse distance d passed by a 
charm quark on the way out of nuclei (in units of its radius), while stars correspond to 
the sum d\ + d 2 of distances for charmed quark and anti-quark. 

d (in units of the nuclei radius R) is shown by a histogram in Fig.^. Note that it is 
basically flat between d~ .2R (or 1 fm for heavy nuclei) and 2R (the diameter). However, 
the distribution of (di + d<i)j R (shown by stars in Fig.|^) is quite different. It is sharply 
peaked at its largest value, but is very strongly suppressed at small ones. In order for both 
charmed quarks to escape, they not only should be created close to the surface, but also 
quarks should be emitted in a very small (tangent) solid angle. As we will show shortly, 
this simple observation is in fact responsible for a significant reduction of the correlated 
charm (and bottom)-induced background. 

A dynamical ingredient of our analysis is dE/dx, the quark energy losses in QGP. 
We would not comment here on a complicated history of its discussion in theoretical and 
phenomenological papers. A consistent treatment (generalizing Landau-Pomeranchuck- 
Migdal approach to QCD) was recently developed by ]7[ (BDPS). The main qualitative 
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difference between QED and QCD cases can brieffy be explained as follows. In QED an 
electron is scattered and has a complicated zigzag-like trajectory, while its field go without 
interaction by a straight line. In QCD it is the quark which is going by approximately 
straight line, while its gluonic field suffer multiple rescatterings. The BDPS result for the 
energy loss is 



where Cr is Casimir for quark color representation, E is the collision energy, \ g is gluonic 
mean free path and /i is the rms momentum transfered in each scatterings. Substituting 
some "reasonable" parameters of QGP at RHIC (corresponding to "hot glue scenario, see 
H H) we have estimated dE/dx ~ 2GeV/ fm. 

Our next step is Monte Carlo simulation of charm production. In order not to intro- 
duce any additional points of discussion, we follow ref.|]IJ] as close as possible. We have 
ignored the "thermal charm" and assumed that each central AuAu collision at RHIC 
produces cc pairs with a very stiff pt distribution^ generated by the leading order and 
resummed next-to-leading-order QCD processes. This approach works in pp case, and 
that is why GMRV has found that the correlated charm decay contributes so strongly at 
large dilepton masses. 

However, after the energy losses dE/dx are included, only very few of c or b quarks can 
in fact escape, while most of them are stopped. Eventually, those should have p t spectra 
similar to all other hadrons, governed by low decoupling temperature T ~ 140MeV and 
hydro effects. Since both thermal and hydro velocities are not large, we have ignored 
them. 

We have simulated semileptonic decays of c and b quarks and show the resulting 
invariant mass M = (pi + +p;_) 2 spectrum in Figs.^, |j. In both cases the histogram shows 
free decays, while stars include the effect of dE/dx. Those two cases are very different: 
while in free space the invariant mass distribution has a smooth and large tail toward 
the large masses, with dE/dx one clearly see two distinct components: charm decay at 
rest and the contribution of escaping ones. The boundary between two components is at 
Mi+i- » 1.7GeV for c and 4.5 GeV for b decays. Above it we have found a background 
suppression, roughly by about 2 orders of magnitude. These features survive reasonable 
modification of charm production spectra or of the chosen dE/dx value. 

How important this reduction may be in practice? In order to answer this question, 
one has to evaluate dilepton production, both primary (known as the Drell-Yan process) 
and secondary (non-equilibrium || and thermal ||) ones. In this short note we will not go 

3 An approximate parameterization used is dN/dpf ~ l/{pt + 0.5) 2 2 , so the power of pt is close to 4. 
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Figure 3: The distribution of dilepton invariant masses produced by a semileptonic 
decays of charmed quarks, with (stars) and without (histogram) the matter effect due to 
dE/dx. 

into discussion of it| and simply return to GMRV estimates. As seen from Fig.[l] the ratio 
(dilepton yield) /(correlated charm background) is about 1/10 for M = 2 — 8GeV, while 
(dilepton yield) /(b decay background) is about 1/3 for M > 5GeV. Those are exactly 
the mass regions where our suppression discussed above appears! Thus we conclude that 
(dilepton yield) /(correlated charm background) is probably above 1; and that b-decays 
are simply negligible. More quantitative conclusion is difficult to get now: also one should 
consider acceptance of the particular detector, etc. 

Since we are still in a situation with the signal/background ratio being around 1, 
additional experimental tools are needed in order to separate dileptons from charm decays. 
At least two are available: (i) dileptons are produced back-to-back in azimuthal angle, 
while leptons from charm decay are nearly isotropic in it; (ii) Drell-Yan pairs have the well 
known (1 + cos 2 9) distribution where 9 is the polar angle between the dilepton direction 
in its CM frame and the beam. Also DY and direct charm should have simple scaling 
^4/3 f rom light nuclei (or peripheral collisions), so any excess over it is an indications for 
secondary processes. 

Let us only mention that GMRV make a very good job on DY, but do not include the non-equilibrium 
one. Also, they treat thermal in the leading order only. Both effects are expected to increase the secondary 
production substantially. 
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Dilepton Mass from B decays 




Figure 4: Same as the previous figure, but for b quark decays. 

In summary: in contrast to GMRV, we think that c and b quarks produced in high 
energy heavy ion collisions should be trapped in matter with very high probability. As 
a result, the background due to correlated semileptonic charm decay does not dominate 
the dilepton spectra for invariant masses above 2 GeV. Optimistically, by using various 
angular distributions, one may probably measure both dileptons and charm. 
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